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ABSTRACT 

We report the detection of interstellar hydrogen deuteride (HD) toward the 
supernova remnant IC443, and the tentative detection of HD toward the Herbig 
Haro objects HH54 and HH7 and the star forming region GGD37 (Cepheus A 
West). Our detections are based upon spectral line mapping observations of 
the R(3) and R(4) rotational lines of HD, at rest wavelengths of 28.502 and 
23.034 /im respectively, obtained using the Infrared Spectrograph onboard the 
Spitzer Space Telescope. The HD R(4) /R(3) line intensity ratio promises to be 
a valuable probe of the gas pressure in regions where it can be observed. The 
derived HD/H 2 abundance ratios are (1.19±g;||) x 10" 5 , (1.80^) x 10~ 5 , and 
(1.4llo; 3 3) x 10~ 5 respectively (68.3% confidence limits, based upon statistical 
errors alone) for IC443 (clump C), HH54, and HH7. If HD is the only significant 
reservoir of gas-phase deuterium in these sources, the inferred HD/H2 ratios are 
all consistent with a gas-phase elemental abundance [nn/^Hlgas ~ 7.5 x 10~ 6 , a 
factor 2-3 below the values obtained previously from observations of atomic 
deuterium in the local bubble and the Galactic halo. However, similarly low 
gas-phase deuterium abundances have been inferred previously for molecular gas 
clouds in the Orion region, and in atomic clouds along sight-lines within the 
Galactic disk to stars more distant than 500 pc from the Sun. 

Subject headings: ISM: Molecules — ISM: Abundances — ISM: Clouds - molec- 
ular processes 
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1. Introduction 

With a cosmic abundance > 10~ 5 , deuterium is more common than all but a half-dozen 
heavy elements. The initial deuterium abundance was set by primordial nucleosynthesis 
(e.g. Schramm & Turner 1998), roughly 100s after the Big Bang, and has been reduced by 
stellar nuclear reactions over the subsequent 13.7 Gyr. The primordial deuterium abundance 
revealed from intergalactic absorption components (e.g. Kirkman et al. 2003, who obtained 
[D/[H]=2.78lg;3g x 10_5 toward Q1243+3047) is in good agreement with models for pri- 
mordial nucleosynthesis, given the baryon densities derived independently from an analysis 
of the cosmic microwave background (e.g. Spergel et al. 2003). Thus, the distribution of 
deuterium in the Galaxy probes both stellar processing and the degree to which material is 
mixed efficiently within the interstellar medium. The abundance of interstellar deuterium 
has been measured primarily by ultraviolet absorption line observations of atomic deuterium 
within diffuse clouds (reviewed by Moos et al. 2002). Recent measurements of the Galactic 
deuterium abundance reveal significant variations (over the range ~ 0.7 — 2.2 x 10~ 5 ) that are 
inconsistent with current predictions from Galactic chemical evolution models (e.g. Wood 
et al. 2004, Friedman et al. 2006, and references therein); understanding the nature of these 
inconsistencies is critical to our understanding of chemical evolution within the Galaxy. 

While deuterium has been extensively studied in atomic clouds, the deuterium abun- 
dance in molecular clouds is less certain. In diffuse molecular clouds, interstellar hydrogen 
deuteride (HD) has been widely observed by means of ultraviolet spectroscopy - starting 
over 30 years ago with the Copernicus satellite (reviewed by Spitzer & Jenkins 1975) - and 
is the most abundant deuterium-bearing molecule. However, deriving the deuterium abun- 
dance from the HD column in diffuse clouds is difficult for several reasons: the chemistry is 
complex, HD contains only a trace amount of D, and the HD abundance is sensitive to the 
density and UV field in the cloud (Lacour et al. 2005; Le Petit et al. 2002). 

In dense molecular clouds, by contrast, the observations are much fewer. Prior to the 
results reported here, detections of hydrogen deuteride emissions had been reported toward 
just one cloud: the Orion Molecular Cloud (Wright et al. 1999; Howat et al. 2002; Bertoldi 
et al. 1999; the latter's detection of HD from warm shocked gas having been anticipated the- 
oretically by Timmermann 1996). In addition, HD R(0) line absorption has been detected 
toward two far-IR continuum sources: Sgr B2 (Polehampton et al. 2002) and W49 (Caux et 
al. 2002). HD is not easily detected by infrared emission-line spectroscopy; its small dipole 
moment leads to relatively weak emission and its large rotational constant places its pure 
rotational transitions at infrared frequencies that are inaccessible from ground-based tele- 
scopes. On the other hand, a large number of other deuterium-bearing molecules are readily 
detected with ground-based radio telescopes. Many molecules, such as methanol, show high 
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levels of deuterium fractionation, in which the deuterated/non-deuterated abundance ratio 
can exceed the cosmic deuterium abundance by more than four orders of magnitude; several 
doubly- and even triply-deuterated species have been observed (Parise et al. 2004). 

In this paper, we report a detection of HD in shocked molecular gas within the supernova 
remnant IC443, clump C, and tentative detections in three other sources - the Herbig-Haro 
objects HH54 and HH7, and the star-forming region GGD37 (Cep A W) - all from observa- 
tions with the Infrared Spectrometer (IRS) on the Spitzer Space Telescope. The observational 
results are discussed in §2, below, and the implied gas density and HD abundance derived 
in §3. The implications of the inferred HD abundance are discussed in §4. 

2. Observations and results 

Two pure rotational lines fall within the 19.5 - 37.2 /im wavelength range covered 
by the Long-High module on the IRS: the HD R(3) and R(4) transitions at 28.502 and 
23.034 /im, respectively 1 . Weak spectral features have been observed serendipitously at 
both wavelengths in several sources, as shown in Figure 1. These plotted spectra were 
obtained in spectral line mapping observations of IC443C conducted in GO cycle 1 - and 
of HH54, HH7 and GGD37 carried out in the IRAC- and IRS-GTO programs - along with 
observations with the Short-High module. In addition, IC443C, HH54, and HH7 had been 
observed in the Short-Low module, to provide complete spectral coverage from 5.2 - 37 fim 
so that the H 2 S(0) - S(7) transitions - and the fine structure lines of [Fell], [SI], [Nell] 
and [Sill] - could all be mapped. These other observations, especially those of the H 2 pure 
rotational transitions, tightly constrain the temperature and column density of the shocked 
gas, allowing a determination of the HD/H 2 abundance ratio. 

Table 1 summarizes the measured R(3) and R(4) line intensities. At each of six positions 
- 1 each in IC443C, HH7 and GGD37, and 3 in HH54 - we have obtained averages of all 
spectra observed within a circular region, the contributing spectra being weighted by a 
Gaussian taper (HPBW of 15") from the center of each such region. The central positions 
are given in Table 1. Full details of the observations toward HH7 and HH54, as well as the 
data reduction methods we have developed, are given by Neufeld et al. (2006). 

Searching the NIST database of atomic fine structure lines - along with the JPL molec- 



1 HD R(0) R(l), and R(2) lie at wavelengths inaccessible to the Long-High module on Spitzer/IRS, while 
R(5) and higher transitions are expected to be so weak that upper limits from Short-High observations fail 
to provide useful constraints. 
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ular line list - for alternative identifications of the 28.502 and 23.034 fim features, we found 
no plausible candidates besides HD. In particular, a large number of water rotational transi- 
tions lie in the wavelength region of interest, but the only such transitions within one spectral 
resolution element of 28.502 and 23.034 /im are very-high-lying transitions that would be 
accompanied by other - much stronger - transitions that are absent in the observed spectra. 
As a check upon the reality of the observed spectral features, we have constructed maps 
of their distribution. Figure 2 shows the map obtained toward IC443C, which exhibits the 
highest column density of warm H 2 and the strongest 28.502 and 23.034 /jm features. Here, 
we compare maps of the 28.502 and 23.034 fxm features with those of the H 2 S(2), which has 
an upper state energy E v /k =1682 K, similar to those of HD R(3) and R(4) (E v /k = 1271 
K and 1895 K). A map of the H 2 S(3) - for which Eu/k = 2504 K - is also shown. The agree- 
ment of the observed morphology in the four lines lends strong support to the identification 
of HD and eliminates the possibility that the 28.502 and 23.034 fim features are (previously 
unidentified) instrumental artifacts. In the other sources, the signal-to-noise ratio is insuffi- 
cient to allow the distribution of the 28.502 and 23.034 /zm features to be mapped reliably; 
accordingly, we conservatively describe the detection of HD in those sources as tentative. 



3. HD abundances obtained with Spitzer 

Using Spitzer /IRS observations of the H 2 S(0) through S(7) transitions towards these 
sources, we can constrain the column density and temperature of the warm, shocked molecu- 
lar hydrogen very well, and can thereby estimate the HD/H 2 abundance ratio. As described 
by Neufeld et al. (2006), we have fitted the H 2 rotational diagrams with a model which 
invokes two components: a warm component at temperature T w ~ 400 K and a hot compo- 
nent at temperature T h ~ 1000 K. The temperatures and column densities of the warm and 
hot gas components are given in Table 1. Our analysis neglects the effects of subthermal 
excitation for the H 2 transitions that we have observed; while possibly important for the 
higher-J transitions of H 2 (i.e. S(6) and S(7)), such effects are negligible for the lower-J H 2 
transitions that probe the most of the gas capable of producing HD R(3) and R(4) emissions. 

Unlike the lower-J states of H 2 , however, the HD level populations for J = 5 and 4 are 
expected to show departures from LTE, due to the presence of a small but non-zero dipole 
moment. We used a statistical equilibrium calculation to compute the HD level populations, 
adopting the collisional rate coefficients of Flower & Roueff (1999) for excitation by H 2 . We 
thereby determined the resultant HD emission, given the two-component model parameters 
derived from our fit to the H 2 rotational diagram, and assuming (1) that HD is at the same 
kinetic temperature as H 2 ; (2) that the HD abundance relative to H 2 , n(HD)/n(H 2 ), is the 



- 5 - 



same in both the warm (~ 400 K) and hot (~ 1000K) gas components; and (3) that the warm 
and hot components have the same gas pressure, p = n(H 2 )T. Typically, the contributions 
made by the hot and warm components to the observed HD emissions were roughly equal 
in our excitation model. The critical densities for R(3) and R(4) - i.e. the densities at 
which the departure coefficients for J = 5 and J = 4 equal one-half - correspond to pressures 
~ 4 x 10 7 and ~ 10 8 cm" 3 K, respectively, so the R(4)/R(3) ratio is a useful density indicator 
for pressures in the ~ 10 7 — 10 8 cm~ 3 K range. The R(4)/R(3) intensity ratios predicted for 
IC443C are 0.32, 0.40, 0.50, 0.61, 0.72 and 0.82 respectively for logi (n[H 2 ]T/cm~ 3 K) = 7.0, 
7.2, 7.4, 7.6, 7.8, and 8.0. 

We varied the gas pressure, p, and the HD abundance relative to H 2 , n(HD)/n(H 2 ), 
to obtain the best fit to the observed HD R(3) and R(4) line strengths. Because the HD 
transitions are optically-thin, the HD line intensities scale linearly with n(HD)/n(H 2 ). While 
our estimates of the gas pressure or density scale in inverse proportion to the assumed 
collisional rate coefficients, the derived HD abundance would be entirely unaffected by a 
(uniform) change in the adopted rate coefficients. Table 1 lists the best-fit gas pressure and 
HD abundance relative to H 2 derived for each region. 

In Figure 3, we show the confidence regions for the gas pressure and HD/H 2 ratio: 
68.3, and 95.4 % confidence regions are shown for each of the five sources we observed, 
based upon the statistical errors on the HD R(3) and R(4) line fluxes and assuming the H 2 
column densities and temperatures (N w , T w , N^, T^) listed in Table 1. The confidence limits 
presented here are based solely upon the statistical errors in the measured line fluxes and do 
not include systematic uncertainties in the HD excitation model. Such uncertainties are hard 
to estimate quantitatively, and include possible errors in the collisional rate coefficients and 
in the assumption that the warm and hot gas components share a common pressure. The 
derived HD abundances (although not the gas pressures) appear to be relatively insensitive 
to the assumed collisional rate coefficients. If the HD-He (Roueff & Zeippen 2000) or HD H 
rate coefficients (Flower & Roueff 1999) are adopted in place of the HD-H 2 rate coefficients, 
the best-fit HD abundances change by less than ~ 10%. On the other hand, if it is assumed 
that the warm and hot gas components share a common density rather than being in pressure 
equilibrium, the best-fit HD abundances increase by up to ~ 50%. 



4. Discussion 

The signal-to-noise ratio is significantly better toward IC443C, HH54FS, and HH7, the 
only cases for which R(4) is detected at the 5 a level. The best fit gas pressures in these 
sources, found to lie in the range 1.6 x 10 7 — 1.0 x 10 8 cm" 3 K, correspond to gas densities 
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in the range 4 x 10 4 to 2.5 x 10 5 cm 3 for the warm gas component. These values are in 
reasonable agreement with those derived previously for these regions (see discussion of HH7 
and HH54 by Neufeld et al. 2006; and of IC443 by Snell et al. 2005.) The HD R(4)/R(3) 
intensity ratio promises to be a valuable probe of the gas pressure in regions where it can be 
observed. 

Our analysis yields best estimates of (1.19+°;^) x 10~ 5 , (1.80+°;^) x 10" 5 , and (1.41+°; 4 ^) x 
10~ 5 (68% confidence intervals) respectively for the HD abundance relative to H 2 in IC443C, 
HH54FS, and HH7, but all three estimates are consistent (within the 68% confidence in- 
tervals) with n(HD)/n(H2) = 1.5 x 10~ 5 . If HD accounts for all the deuterium in the 
gas phase 2 , this would imply [nD/nn]gas = 7.5 x 10~ 6 , a value which lies a factor 2 - 
3 below those inferred from atomic D absorption line observations of the local bubble 
([n D /n H ] g as = (1.52±0.08) xlO -5 ; Moos et al. 2002), atomic D absorption line observations of 
the halo ([n D /nu] ga s = (2.2 ±0.7) x 10~ 5 ; Sembach et al. 2004), and the expected abundance 
from primordial nucleosynthesis ([ud/^h] = (2-62+q;2q) x 10~ 5 ; Spergel et al. 2003, given the 
cosmological parameters determined by the Wilkinson Microwave Anisotropy Probe). How- 
ever, it is comparable to those derived for the Orion Molecular Cloud from a claimed detection 
of HD R(5) toward Orion Peak 1 ([n D /n H ]ga S = (7.6 ± 2.9) x 10~ 6 ; Bertoldi et al. 1999) and 
from observations of HD R(0) emission from the Orion Bar {[n^/n-[i\ gas = (1.0 ±0.3) x 10~ 5 ; 
Wright et al. 1999). It is also remarkably consistent with the deuterium abundances observed 
recently in atomic clouds along long sight-lines within the Galactic disk toward stars further 
than 500 pc from the Sun and with atomic H column densities in excess of 10 20 5 cm -2 , which 
yield [n D /n H ] g as = (8.5 ± 1.0) x 10~ 6 (Hoopes et al. 2003; Wood et al. 2004). 

The explanation of these variations has been a subject of considerable debate (Friedman 
et al. 2006). Hebrard and Moos (2003) have argued that the lower values inferred for high-A^ H 
sight-lines to distant stars represent the typical Galactic elemental abundance for deuterium, 
requiring a substantial degree of destruction via astration (i.e. in stars), while Draine (2004) 
has argued for the importance of depletion in reducing the gas-phase deuterium abundance, 
a possibility originally raised by Jura (1982). In this picture, material within the local 
bubble is anomalous because it has suffered recent grain destruction in shocks driven by 
stellar winds or supernovae. An observed correlation between the gas-phase abundances of 



2 In interpreting their observations of Orion, Bertoldi et al. argued that significant destruction of HD 
occurred in the shocked gas that they observed, as the result of the reaction H + HD — > D + H2 (Timmermann 
1996), and applied a correction in deriving the gas-phase elemental abundance of HD. In the shocked regions 
that we have observed, however, the inferred shock velocities (~ 10 — 20kms~ 1 ; Neufeld et al. 2006) are 
significantly smaller than those for the Orion shock; thus the resultant atomic H abundance is expected to 
be much smaller and the extent of HD destruction negligible. 
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deuterium and the refractory element titanium (Prochaska et al. 2005), and the possible 
detection of deuterated polycyclic aromatic hydrocarbons (Peeters et al. 2004), have been 
interpreted as supporting the importance of depletion. While our determination of the HD 
abundance in dense shocked clouds does not directly settle this debate, it does suggest that 
the gas-phase deuterium abundance in dense clouds differs little from that inferred for typical 
diffuse clouds along high-A^n sight-lines to distant stars within the Galactic disk. 

This work, which was supported in part by JPL contract 960803 to the Spitzer IRS 
Instrument Team and by RSA agreement 1263841, is based on observations made with the 
Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California 
Institute of Technology, under a NASA contract. D.A.N, and P.S. gratefully acknowledge 
the additional support of grant NAG5-13114 from NASA's LTSA Research Program. 
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Table 1. Observations of HD 
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19.84 
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T h (K)» 


N/A 


1047 


1029 


1150 


1062 


1025 


logio (Affe/cm- 2 )" 


N/A 


20.01 


19.17 


18.86 


19.13 


19.22 


R(3) intensity / 10 -7 ergcm~ 2 s -1 sr -1 


29.0±2.8 C 


52.3 ± 3.1 


17.9 ±2.2 


8.0 ±1.2 


15.6 ±2.9 


17.7 ± 1.0 


R(4) intensity / 10~ 7 crgcm~ 2 s -1 sr -1 


13.6 ±2.6 


21.7 ±3.1 


8.7 ± 1.2 


6.2 ± 1.3 


5.8 ± 1.7 


9.3 ± 1.8 


7(R(4))//(R(3)) 


0.469 ±0.100 


0.415 ± 0.064 


0.486 ±0.092 


0.770 ± 0.199 


0.374 ±0.130 


0.527 ±0.107 


log 10 (n(H 2 )T/cm- 3 K) (best fit) 


N/A 


7.23 


7.55 


7.98 


7.22 


7.53 


log 10 (n(H 2 )T/cm- 3 K) (68.3% c.l.) 


N/A 


7.10 - 7.36 


7.37 - 7.75 


7.64 - 8.71 


6.88 - 7.51 


7.32 - 7.73 


log 10 (n(H 2 )T/cm- 3 K) (95.4% c.l.) 


N/A 


6.93 - 7.49 


7.20 - 8.03 


> 7.33 


6.53 - 7.89 


7.06 - 7.98 


10 s x ji(HD)/n(H 2 ) (best fit) 


N/A 


1.19 


1.41 


0.53 


2.07 


1.80 


10 s x ji(HD)/n(H 2 ) (68.3% c.l.) 


N/A 


0.95 - 1.54 


1.07 - 1.87 


0.41 - 0.71 


1.20 - 4.46 


1.49 - 2.34 


10 5 x n(HD)/n(H 2 ) (95.4% c.l.) 


N/A 


0.79 - 2.23 


0.82 - 2.68 


0.33 - 1.02 


0.72 - 9.93 


1.27 - 3.62 



a Two-component fit to the H 2 line intensities with a warm component at temperature T w and H 2 column density N w , and a hot component at 
temperature T), and H 2 column density 

b not available, because the H 2 S(3) — S(7) lines have not been observed 

c la statistical uncertainty 
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Fig. 1. — HD R(3) and R(4) spectra observed for 15" (HPBW) diameter circular apertures 
centered at 6 locations defined in Table 1. From top to bottom the spectra are those of 
GGD37, IC443C, HH7, HH54C, HH54E+K, and HH54FS. Arbitrary continuum offsets have 
been introduced for clarity. The strong feature at 28.2188 jum is the S(0) line of H2. 




Fig. 2.— HD R(3), HD R(4), H 2 S(2) and H 2 S(3) line intensities observed toward IC443C. 
Logarithmic contours are shown, in steps of 0.1 dex, with the lowest contours corresponding 
to line intensities of 2 x 10~ 6 , 1 x 10~ 6 , 2 x 10~ 4 and 1 x 10~ 3 ergcm~ 2 s _1 sr _1 respectively 
for the HD R(3), HD R(4), H 2 S(2) and H 2 S(3) transitions. The rectangular boxes demark 
the mapped regions for the Short-Low (H 2 lines) and Long-High (HD lines) modules. The 
effective spatial resolution is ~ 3" and ~ 10" (HPBW) respectively for the H 2 and HD 
transitions. The horizontal and vertical axes show the R.A. (Aacos5) and declination 
(AS) offsets in arcsec relative to a =06h 17m 44.20s, 5 = 22d21'49.1" (J2000). The spectra 
presented in Figure 1 are averages over a 15" region centered at offset (—24", —20"). 
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Fig. 3. — Contours for A% 2 = 1 (inner contour) and 4 (outer contour) as a function of 
gas pressure, n(H 2 )T, and the HD abundance relative to H 2 , n(HD)/n(H 2 ), are shown for 
each of the five sources we observed. These results adopt the H 2 column densities and 
temperatures (N w , T w , Nh, Th) listed in Table 1 as priors, and assume that the observational 
errors in the measured line intensities have a Gaussian distribution function. The projections 
of these contours onto the horizontal and vertical axes yield the 68.3% (inner contour) and 
95.4% (outer contour) confidence limits on n(H 2 )T and n(HD)/n(H 2 ). The joint probability 
contents for the inner and outer contours are 39.3 and 86.5 % respectively. Note the different 
range of gas pressures and HD abundance shown in each panel. Crosses denote the best-fit 
values (Ax 2 = 0). 



